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Abstract 

We establish the convergence of statistically invariant states of the three dimensional Boussinesq 
Equations in the infinite Prandtl number limit. The equations are subject to a temperature gradient on 
the boundary and to internal heating in the bulk driven by a stochastic, white in time, gaussian forcing. 

For the active scalar equations given by the infinite Prandtl system we show that the associated invariant 
measures are unique and that moreover the corresponding Markovian dynamics evolving probability dis¬ 
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1 Introduction 

Buoyancy driven convection plays a central role in a wide variety of physical processes: from earth’s cli¬ 
mate system to the internal dynamics of stars. As such it is of fundamental importance to identify and 
predict robust statistical quantities in these complex flows and to connect such statistics with the basic 
equations governing their dynamics, for example the Boussinesq equations. In particular characterizing pat¬ 
tern formation, mean heat transport, and small scale dynamics as a function of physical parameters and 
boundary conditions remains a topic of intensive research theoretically, numerically, and experimentally; see 
e.g. [BPAOO, Man06, AGL09, LXIO] for a broad overview of recent developments. 

It has long been understood that statistically invariant states of the nonlinear partial differential equa¬ 
tions of fluid dynamics provide mathematical objects which are expected to contain various robust statistical 
quantities found in turbulent fluid flows. An ongoing challenge is therefore to address the existence, unique¬ 
ness, ergodicity, and dependence of these measures on parameters in a variety of specific contexts. While 
one may certainly pose these questions for deterministic equations (cf. [FMRTOl]) the stochastic setting can 
be more tractable given the regularizing effect of noise on the associated probability distribution functions. 
Moreover stochastic terms can be physically reasonable in various settings. In convection problems energy 
may be supplied to the system through both boundary and volumetric heat fluxes, which for instance models 
radioactive decay processes in the earths mantle; see [Rob67, TZ67, LDB04, WDll, GS12, BN12]. Both 
sources can therefore have an essentially stochastic character in situations of physical interest. 

In this and a companion work, [FGHR] , we study statistically invariant states of the stochastically driven 
Boussinesq equations 


-f u • Vu) — Au = Vp -I- i?akr, V • u = 0, (1.1) 

Fr 

N 

dT FvL-yTdt = + (1.2) 

fe=i 

for the (non-dimensionalized) velocity field u = (ui, 142 ,^ 3 ), pressure p, and temperature T of a buoyancy 
driven fluid. The system (1.1)-(1.2) evolves in a three dimensional domain V = [0,L]^ x [0,1] and is 
supplemented with the boundary conditions 

U| 2 ^q = = 0, T\z=o = T\z=i =0; u,T are periodic in x = (a:i,a; 2 ). (1-3) 

The dimensionless physical parameters in the problem are the Prandtl number Pr and Rayleigh numbers Ra 
and Ra] see Remark 3.3 below for further details. The driving noise is given by a collection of independent 
white noise processes dW^ = dW^ (t) acting in spatial directions Cfc = at {x) which form a complete orthogo¬ 
nal basis of eigenfunctions (ordered with respect to eigenvalues) of the Laplace operator on T) supplemented 
with homogeneous Dirichlet boundary conditions for z = 0,1 and periodic in x. The stochastic terms in 
(1.2) have been normalized so that 


N 

ll<^fe|lL2(r>) = 1J 

fc=l 

with the strength of the body forcing expressed in the physical parameters Ra and Ra] see (3.7) below. 

Our principal aim here will be to establish convergence properties of statistically invariant states of (1-1)- 
(1.3) in the Large Prandtl number limit when the physical parameter Pr in (1.1) diverges to 00 . In other 
words we will establish the convergence to invariant measures of the active scalar equation 

— Au = Vp + RakT, V • u = 0, 

N 

dT + u- VTdt = ATdt + akdW'" 

k=l 

complemented with boundary conditions for T and u as in (1.3). The analysis of convection in the large 
Prandtl number limit is of basic interest in a variety of physical contexts, most notably in modeling certain 


(1.4) 

(1.5) 
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portions of the earth’s mantle and for convection in gasses under high pressure, where the Prandtl number 
can reach the order of 10^^, see [CD99, DCOl, OSll]. It is worth emphasizing that the system (1.4)-(1.5) 
has very complex dynamics even without stochastic forcing when the Rayleigh number(s) are sufficiently 
large; see [CD99, BPAOO, DCOl, Wan04a, Par06, AGL09, LXIO, OSll], 

Let us now present a heuristic version of our main results to provide intuition. For the precise formulation 
see Theorem 3.1 below. 

Theorem 1.1. Fix any Ra,Ra > 0 and consider (1.1)-(1.3) and (1.4)-(1.5) with N > N{Ra, Ra) indepen¬ 
dently forced directions in the temperature equation. Then (1.4)-(1.5) possesses a unique ergodic invariant 
measure poc- Let {ppr}pre(o,oo) be a sequence of statistically invariant states associated to (1.1)-(1.3).^ 
Then ppr converges to p^o suitable metric. In particular, for any sufficiently regular observable <f) on 
the phase space of T, 


4>{T)dppr - j 4>{T)dpa 


< C{Pr)-<^ , 


( 1 . 6 ) 


where C = C[4), Ra, Ra), q = q{Ra, Ra) > 0 are independent of Pr and q is independent of 4>. 

Our method of proof can be simplified to yield analogous results for the 2D version of the stochastically 
forced Boussinesq equation (1.1)-(1.3). Theorem 1.1 may also be seen as a complement to a series of recent 
works [Wan04a, Wan04b, Wan05, Wan07, WanOSa, WanOSb] which addresses large parameter limits for flows 
driven only through the boundary and in a deterministic framework. Here we show that the addition of 
stochastic terms to the Boussinesq system allows for stronger convergence results, but the proofs require a 
different framework. 

The starting point of our analysis is to establish a strict contraction for the semigroup {P°}t>o evolving 
probability distributions associated to (1.4)-(1.5). We will show that there exists k € (0,1) such that for 
> 0 sufficiently large 


Pil^Pt ,: P) (1-7) 

for any measures p, p on this phase space associated with T. Here p is a Kantorovich-Wasserstein metric 
defined by a carefully chosen distance on the phase space; see (2.9) below. In a series of recent works, 
[HM08, HMIO, HMSll], a general scheme for proving a strict contraction, as in (1.7), was developed. 

Once an estimate of the form (1.7) is established, the convergence of statistically invariant states, as in 

(1.6) , can be translated to the convergence of solutions on finite time intervals. More specifically, since p 
is a metric defined by Lipschitz norms of observables, when vprft), T^'^{t) and Voo{t), T°°{t) represent the 
probability laws and associated solutions of (1.1)-(1.2) and (1.4)-(1.5) respectively, one may show that 

p{>^Pr{t),Poo{t)) < C (Eexp(p||r^’'(t)f). Eexp(p||T-(t)f) • E||T^’'(t) - T^{twf'' (1.8) 
for suitable choices ol C,r],q > 0 independent of Pr,t > 0. By then combining (1.7) with (1.8) we find that 

p(n/ip„po)<c(Eexp(p||rJ’'(t,)f)-Eexp(p||T“’^’-(t,)f).E||Tj’'(t,)-r|°’^"(i*)ir)'^" (1-9) 

where T^’’ are stationary solutions of (1.1)-(1.2) and T^P^ are solutions of the infinite Prandtl number 
system starting from initial conditions given by Tg'^{0). Here H is the projection onto the T portion of 
the phase space and the time t* in (1.9) comes from (1.7). For specific details on the relationship between 

(1.7) -(1.9) see Section 2.2 below. 

In summary through (1.9) the strict contraction property (1.7) reduces the question of convergence 
of stationary states to establishing some Pr-independent exponential moment bounds, and proving the 
convergence of solutions of the finite Prandtl system to those of the infinite Prandtl system at a fixed positive 
finite time. Moreover (1.9) provides a means for establishing the rate of convergence for p{Ilppr, po) —?► 0 

^In the precise statement of this result, Theorem 3.1, we will impose an additional technical assumption on ppr- This 
condition is roughly analogous to a finite energy criteria for weak solutions of the 3D Navier-Stokes equations. In [FGHR] we 
have demonstrated the existence of such states UPr- 
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as Pr oo. It is also worth noting here that the approach leading to the bound (1.9) does not rely on 
uniqueness of statistically invariant states for finite values of Pr >0.^ 

We therefore face two challenges. The first is to establish the strict contraction (1.7) for (1.4)-(1.5). 
Guided by the classical Doob-Khasminskii Theorem [Doo48, Km60, DPZ96] and as encompassed by the 
more recent developments in [HMDS, HMIO, HMSll] one would expect a contraction of the type (1.7) when 
the Markov semigroup is smoothing, suitable moment bounds hold, and there is some form of irreducibility 
in the dynamics. The second major challenge is to address to the finite-time convergence desired in (1.9). 
We stress here that this convergence to zero does not follow from continuous dependence on parameters. 
Indeed, the system (1.1)-(1.3) with large Pr is a singular perturbation of (1.4)-(1.5). For the infinite Prandtl 
system (1.4)-(1.5), u is determined from T at each time step by solving a Stokes problem and thus u does 
not satisfy an independent evolution equation. On the other hand one prescribes in (1.1)-(1.3) an initial 
condition for u which is absent in the limit Pr — >■ oo. 

Regarding the hrst challenge, the question of smoothing for the Markov semigroup can be translated to 
a control problem; see (B.IO) below. In our setting, when the number of forced directions N = N{Ra, Ra) 
is sufficiently large, an appropriate control can be found through Foias-Prodi type considerations [FP67]. 
Since (1.4)-(1.5) may be seen as an advection diffusion system with u being two derivatives smoother than 
T, such a strategy largely repeats the approach used in previous works on the 2D stochastic Navier-Stokes 
equations [EMSOl, HM06, HMDS, KS12]. On the other hand establishing suitable moment bounds is more 
delicate due to the non-homogenous boundary conditions imposed in (1.3) and requires a careful use of the 
maximum principle along with exponential martingale estimates. These bounds have been carried out in the 
companion work [FGHR] . 

Thus the main obstacle to demonstrating (1.7) is to establish irreducibility, which is not so transparent 
from the path set out in previous works e.g. [EMOl, HM06, GGIIV13, FGHRT13]. This is because the 
system (1.4)-(1.5) with its stochastic terms removed can have highly non-trivial dynamics; [CD99, DCOl, 
Wan04a, Par06, OSll]. We proceed to show that despite this observation every invariant measure of (1.4)- 
(1.5) contains zero as a common point of support. Indeed we establish with the use of another Foias-Prodi 
bound that a Girsanov shift of (1.4)-(1.5) has almost trivial dynamics with positive probability. We then 
employ moment estimates and stopping time arguments to translate this non-zero probability back to the 
original system (1.4)-(1.5) yielding the desired irreducibility 

We turn to the second major challenge: the convergence of solutions at a fixed finite time as Pr —?► oo, 
as desired in (1.9). Here we need to develop a suitable asymptotic analysis for (1.1)-(1.2). This is a 
non-trivial task since the small parameter lies in front of the time derivative terms in (1.1). Similarly to 
[Wan04a, Wan04b, Wan05, Wan07, WanOSa, WanOSb] we derive an ‘intermediate system’ which is close 
to the finite Prandtl system over bound time intervals and which converges to the infinite Prandtl system 
after an 0{e) time transient. However, in our case time regularity properties crucially used in previous 
works is missing due to the stochastic forcing. As a substitute we derive a stochastic evolution equation 
for the velocity component and use martingale properties of associated Ito integrals. Our analysis then 
takes advantage of some previously unobserved cancellations in the error terms and delicate stopping time 
arguments. 

The manuscript is organized as follows. We begin with a ‘warm-up problem’ by carrying out an asymptotic 
analysis for a system of stochastic ordinary differential equations which model some of the difficulties which 
are encountered in the analysis of the full system. This preliminary step gives us the occasion to introduce 
the formalities of the Kantorovich-Wasserstein metric in Section 2.2. In Section 3 we introduce the rigorous 
mathematical setting of Boussinesq equations, (1.1)-(1.3), which will serve as a foundation for the rest of the 
analysis. Section 4 is devoted to establishing the contraction (1.7) for the infinite Prandtl system (1.4)-(1.5). 
In Section 5 we carry out the finite time convergence analysis which follows ideas from the formal asymptotic 
procedure derived for the toy model. An appendix recalls various elements essentially contained in previous 
works that we have used in our analysis. Section A is devoted to recalling various moment estimates from 
[FGHR] for a class of drift-diffusion equations and for (1.4)-(1.5). In Section B we detail gradient estimates 
on the Markov semigroup given by (1.4)-(1.5) which are carried out in a similar fashion to e.g. [HM06]. 


^In fact we do not even require (1.1)—(1.3) to define a Markovian semigroup, since we only use the existence of statistically 
invariant states as in [FG95] which are established in [FGHR] via Galerkin truncation in the velocity equation in (1.1)-(1.2). 
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2 A Simple Stochastic Model for the Large Prandtl Number Asymp¬ 
totic 

To understand some of the main issues involved in the large Prandtl number limit we begin by considering 
the simplified setting of a class finite dimensional toy models. Both the usage of contraction properties for 
the Markov semigroup defined by the limit system and the formal asymptotics employed to derive finite time 
convergence results will guide our analysis of the full system (1.1)-(1.3) in Sections 3-5 below. The reader 
anxious to see the analysis of the Boussinesq equations can skip immediately to Section 3. 

2.1 The Equations and General Problem Set-up 

Fix Ml, M 2 , N > 0 and let W = (Wi ,..., Wn) be a standard iV-dimensional Brownian motion, relative to 
some ambient filtered probability space ( 12 , {J^t}i>o,P)- For each e > 0 , consider the following coupled 

systems of stochastic ordinary differential equations 

+ ai{u‘^) = Ra-e{9‘^), 

N 

d9^ + {b2{u^, 0") + 02(0")) dt = Y^ akdW\ 

evolving on Here : fl x [0, 00 ) —>■ 0® : 12 x [0, 00 ) —>■ are adapted stochastic processes 

and (Tfc are fixed elements in . We denote |ctP = fFat is | • | the Hilbert-Schmidt norm on 

N X Ml matrices. The operators —>■ i = 1,2 and e : —>■ are linear and we suppose 

that oi, 02 are strictly positive definite with 

(oi(u),m) > |wp, (O2(0),0) > |0p, (2.3) 

for all u G R-^b 9 G R-^^^ assume that 61 : R^^i —>• R^b 62 : R^i +^2 ]^M 2 g^j.g bilinear and satisfy 

the cancellation properties 

(b 2 (u,u),u) = 0, (b 2 (u, 9 ), 9 ) = 0, (2.4) 

for any u,u G R^^, 9 G R^^. We denote the operator norm of a multilinear operator Q hy \Q\ and note that 
although (2.3) imply |a“^| < 1, we indicate below the dependence of constants on |a“^|. For any fixed value 
of the parameter Ra > 0 we our goal is to study the behavior of statistically invariant states 27| = ('Og,0g) 
of ( 2 . 1 )-( 2 . 2 ) in the limit as e —>■ 0 and its relationship to solutions of 

ai(u°) = Ra-e{9°), 

N 

d9°+ {b 2 {u°, 9°) + 02 ( 0 °)) dt = Y^ (JkdW'^, 0°(O) = 00, 

fc=i 

the formal limit system. 

Under the given assumptions it is not hard to see that, for any fixed values of £ > 0, Ra > 0 the 
problems (2.1)-(2.2) and (2.5)-(2.6), complemented with the initial conditions Uq = (mo,0o) G R^i +'''^2 ^nd 
9q G R^2 respectively, posses unique, probabilistically strong solutions depending continuously on the initial 
conditions; see e.g. [AFS08]. Having fixed Ra > 0, we denote {P(^}t>o, {pn t>o as the associated (Feller) 
Markov semigroups; see (3.11), (3.12) below. Denote Hg : R^i+^^s as the projection onto the 0 

component of {u,9), that is, llg{u,9) = 9. 

The existence of invariant measures for {Pt"}t>o and {P°}t>o easily follows from a standard Krylov- 
Bogoliubov procedure [KB37] and energy type estimates for (2.1)-(2.2) and (2.5)-(2.6). These, possibly 
non-unique, invariant measures are denoted pe and po respectively. If I/q g = (uq g, 0 q g) is an Pb-nieasurable 
random variable distributed as pe, then we denote the correspond solution of ( 2 . 1 )-( 2 . 2 ) as (m|, 0 |) and call 
this a statistically steady state of the system. Denote by Hgpe the push forward of pe under Hg.^ 

^Recall that the push forward is defined as 

UokiAA) = for any A e 


(2.5) 

( 2 . 6 ) 


m "( 0 ) = uo. 

( 2 . 1 ) 

0"(O)=0o, 

( 2 . 2 ) 


du^ , , E E\ 

— -h bi{u ,u ) 


(2.7) 
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2.2 Contraction in the Wasserstein Metric and Its Consequences 

With our toy model now in hand we next recall the general setting of the Kantorovich-Wasserstein distance 
in which we will establish the convergence of statistically invariant states as e —>■ 0. We then give the precise 
details of the strategy outlined in (1.7)-(1.9) in the introduction which, under the additional assumptions 
on (2.5)-(2.6), reduces flg/r® —>• /r° to the question of the closeness of solutions of (2.1)-(2.2) to (2.5)-(2.6) 
on finite time intervals. 

We make use of the Kantorovich-Wasserstein distance as follows: let (X, || • ||) be a Banach space and for 
77 > 0 we consider a weighted metric on X as 

p^{e,e)= inf [ exp(r 7 || 7 f )|| 7 '(s)||ds, for any 6», 0 £ X. 

7eCh[0.1];X)_Jo 

7 (o)=e, 7 (i)=e 


Notice that 


\\0 - e\\ < 0) < eMMPf + l|0f ))l|0 - o\\. (2.8) 

Let Pri(X) be the set of Borel probability measures p onX with f prj{0, U)dp{U) < 00 . On Pri(X) we define 
the Kantorovich-Wasserstein metric, relative to pr^ equivalently as 


Pvih,h) ■= 


sup 

II Lip,r]'^ 1 




inf 


Pr,i0,0)dr{0,0) 


(2.9) 


where ||(/)||Lip,,, := sup^^^for ( 7 !) : X —>• R. The set C{p,fl) is the collection of couplings that is 
Borel probability measures P in Pr(X x X) with p,, p, as its marginals. This mean that the last term in (2.9) 
is equivalent to 


Pr,{p,p) = infEp^(X,r), 


( 2 . 10 ) 


where the infimum is taken over all X-valued random variables X, Y distributed as p^ p respectively. See e.g. 
[Vil08] or [Dud02] for a more detailed treatment and further background on these metrics. 

For any 77 > 0 we define the set of observables 


11^(0)11 +||V<^(0)|| ^ 

14)(X) := <(p € C (X) : 11011^ := sup- . - < 00 


exp(77||0P) 

It is shown in [HM08, Proposition 4.1] that for any (p £ C^(X), 

CUhr, < mL^p,r, < cm,. 


( 2 . 11 ) 


These bounds are useful for translating the convergence of measures in p, to convergence of observables; see 
(2.14) below. 

For simplicity we will now take as an assumption that the semigroup {Pt}t>o associated to (2.5)-(2.6) 
is contractive in p, for some 77 > 0 depending only on fixed constants in the problem, Ra, jcrj. Specifically 
we assume that there exist > 0 and k £ (0,1) such that 


PpipPt., t — ^Prjipj p)} 


( 2 . 12 ) 


for any p,p G Pri(R“2). 

Remark 2.1. In view of results from [HM08, HMIO, HMSll], the condition (2.12) can be verified in many 
interesting situations in finite and infinite dimensions. In particular one may establish this condition (2.12) 
for our toy (2.5)-(2.6) if the Hormander bracket condition is satisfied. On the other hand we verify that 
(2.12) for the Boussinesq system with Pr = 00 , (1.4)-(1.5) in Section 4 below. 
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Under (2.12) the invariant measure po for (2.5)-(2.6) is clearly unique and hence ergodic. Observe that 
for any invariant measure pe of {Pf }t>o, e > 0 , we have that 


P-niJdpei Po) — Pr;(n(/ieP(^ ), pgPf^) < Pjj{Il{pgPf^), (Ilpg)Pf^ ) + Prj{{Tlpe)Pt^ ^ PoPt,) 
< Pr^{IV{pePl), {npe)PtJ + Kpr^iUp,, po) . 


Let 0|(U) = (mq g, 0Q g)) and 0g®(U) = 6*°(U,0o,s) denote respectively solutions of (2.1)-(2.2) and 

(2.5)-(2.6) at time U with initial conditions distributed as p^ and TlePe- Clearly the joint distribution of 
0 g(t^,), dg^(U) has marginals Tl^pgPfJ, (Jlpg)Pf^. 

By (2.10) and (2.8) we have that for any ri,q> 0 


pMPePl),(nPe)Pl) < Ep^i0UU),0l’^{U)) < E (exp(277(|0|(U)P + |0s’"(U)n)|0^(U) - e°/{U)\) 

< CE (exp(3r;(|d|(U)p + |d°’^(U)n)|0|(U) - 0^’^(U)|’/^) 

< c(Eexp(1277|0|(U)n .Eexp(12?7|0°’"(U)n)y^^ (e\0UU) - 0°/ 


where C > 0 depends only on q and p. Combining these two observations and rearranging we infer for any 
p,q>0 

P^(n/ie,Mo)<c(Eexp(12r7|0|(U)n-Eexp(12r;|0°’"(U)n)y^^(E|d|(U)-0°’"(U)l’)'^', (2.13) 

where C depends only on q, p and k. Thus to prove (2.12) it suffices to establishing uniform in £ > 0 moment 
bounds and mean convergence properties on finite time intervals. 

In summary working from (2.13) we now prove. 

Proposition 2.1. Fix any Ra > 0 and consider the systems (2.1)-(2.2) and (2.5)-(2.6) satisfying (2.3), 
(2.4). Assume furthermore that (2.5)-(2.6) satisfies the contraction condition ( 2 . 12 ) for some fixed t^, > 0, 
K G (0,1), and p < g|^| 2 . Let pQ be the unique invariant measure associated to (2.5)-(2.6). Then, for any 
collection {tie}£>o of invariant measures corresponding to ( 2 . 1 )-( 2 . 2 ) one has 


P^{PePe,ho)<Ce'i, (2.14) 

where C = C(\u\, Ra,ttf, k) and q = q{\a\, Ra,t^:, k) > 0 are both independent of the (possibly non-unique) 
choice of elements {/ie}e>o. This implies that for any collection of statically stationary solutions (ug^^s) 
(2.1)-(2.2) and the unique stationary solution 0g of (2.5)-(2.6), 

|E(0(d|)-0(0^))|<C'||<^||2„£'^ (2.15) 

for any observable f) G V 2 r;(E^^) and C = C{Ra, |(t|,U,k) independent of (f and £ > 0. 


2.3 Uniform Moment Bonnds 

The first step in the proof of Proposition 2.1 is to establish uniform moment bounds for solutions of (2.2) 
independently of £ > 0. Such bounds which essentially rely on energy estimates and exponential martingale 
bounds should be expected in view of (2.3)-(2.4). Specifically we prove 

Lemma 2.1. Fix any Ft-adapted process v G C([0, oo);R-^^) and S Fq o,nd consider any adapted f G 
(^([O, oo); R'^i) solving 


N 

df + {b 2 {v, 0 + a 2 {f))dt = Y, ^kdW\ ^( 0 ) = Co , 


(2.16) 


where 02 and 62 are subject to (2.3), (2.4). Then: 
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(i) For any 77 < 2 i^> “*^2^ e,K > 0 and any solution ^ of (2.16) one has 

^ ( t>o > K < e”''^ . 

Consequently, for any t>t), 

Eexp[? 7 [ sup |C(s)P+ [ IC(s)P'^s| I < C'exp(77i|crp)Eexp(?7|5on: 

V \selo,t] Jo J J 

where C = C{r]) is independent of v. 

(a) For any t>0 and rj < 

p ( \m\^ - kl" - e-‘|eoP > K\ J-o) < e-^^, 

and hence for every t > 0 

Eexp (77|^(i)P) < Cexp(77|(7p)Eexp(?7e“‘|^on , 

where C = Cirf) is independent oft and v. 

(Hi) For any stationary solution fs of (2.16) and any tj < 

Eexp(77|^s|^) < C < 00 , 


(2.17) 


(2.18) 


(2.19) 


( 2 . 20 ) 


( 2 . 21 ) 


where C = Cirf) is independent of v. 

(iv) Finally if {Og^^s) ® stationary solution of (2.1)-(2.2) then for any fj < 

Eexp(^|M 5 p) < C < 00 . ( 2 . 22 ) 


where C = C(f)) is independent of £ > 0. 

Proof. By applying Ito’s lemma to (2.16) and using (2.3), (2.4) we obtain d|^p + 2|^p(it < \a\'^dt + 2{^, a)dW. 
Hence with the use of exponential martingale bounds and standard estimates we establish (2.17)-(2.21) as 
in e.g. [HM06, HM08, KS12, GH14]; see also Appendix A below. 

It remains to prove (2.22). To estimate {ug,0g) we multiply (2.1) with Ug. Due to (2.3) and (2.4) we 
have < i?a^|ep| 0 |p implying that 


|u|(t)P< exp(-e ^t)|M|op+ f exp(£ ^(s - t))| 6 »|pds < exp(-e 4)|u| + i?a 2 |ep sup | 6 »|(s)p. 

£ Jo sG[0,t] 


On the other hand 9g satisfies a system of the form (2.16). Thus by (2.18) and (2.21) for any ry < 4 ^^ 2 |e| 2 |q .|2 
and any M > 0 


Eexp( 77 |ug(t)|^ A M) < C'exp(t) (Eexp (27)(exp(—£ ^t)\ugQ\‘^ A , (2.23) 

where C = C{fj) is independent of e > 0 and M > 0. Since the distribution of Ug{t) is independent of t, 
one obtains for sufficiently large t such that exp(—£“^t) < ^ that Eexp( 77 |u|p A M) < C for a constant 
C independent of M > 0,e > 0. The bound (2.22) now follows from the monotone convergence theorem, 
completing the proof. □ 
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2.4 Convergence on Finite Time Intervals 

With the e-independent bounds (2.18)-(2.22) in hand, we next establish convergence on finite time intervals 
to obtain the desired decay in (2.13). Throughout this section for any e > 0 we denote 0g and Ug^, 0g^ 
respectively as the solutions of (2.1)-(2.2) and (2.5)-(2.6) supplemented with initial conditions (uq g,0Q g) 
and 0 g g distributed as a pe and Rpe where pe are the invariant measures of ( 2 . 1 )-( 2 . 2 ). 

Remark 2.2. For clarity of presentation we have restricted our analysis to the initial conditions for (2.1)~ 
(2.2) and (2.5)-(2.6) suitable for establishing (2.13). However, straightforward modifications of forthcoming 
methods show that for any t > 0 and any collection {(moj ^o)}£>o C {itQ}e>o bounded and 

9g —>■ 9g we have 

E|r (t, (ul, 9-g)) - 9\t, 0o)r < C{e'i + % - ^o^), 

from some suitable C = C{Ra,t, |CT|,supg>o I'^ol) > ^ and q = q{Ra,t, |tT|) > 0, both independent of £ > 0. 

2.4.1 Formal Asymptotics: The Corrector 

Our first step is to derive an intermediate system of equations that is close to both the systems (2.1)-(2.2) 
and (2.5)-(2.6) but only after an 0(e) time transient for the limit system. 

To this end we observe that (2.1)-(2.2) is formally approximated by 

d ^ 

—iF + £~^ai(iF) = £~^Ra ■ e(F), dF + {h 2 (u‘^, + a 2 (^^)^ dt = ^ a^dW^. 

fe=i 

Thus, under the (as yet unjustified) supposition that there is a clear separation of time scales between the 
motion of and 9^ we therefore propose the effective dynamics for ( 2 . 1 )-( 2 . 2 ) as 

iF(t) = exp(— e“^ait)MQ s -I- af^{Ra ■ e(9^)) — eyip(—£~^ait)af^(Ra-e{9gg)), (2.24) 

N 

dF+{b2(u^,F) + a2(e^))dt = Y,^kdW'^, F(9)=9l^g. (2.25) 

fc=i 

In simple terms we are making use of the fact that solutions of systems of the form ^u+£~^ Au = £~^g, u(0) = 
Mo, where A is a matrix and 5 is a constant, admits the Duhamel representation u(t) = exp(—£:“^At)Mo + 
(A~^g — exp{—e~^At)A~^g). 


2.4.2 Convergence of the Corrector to the Limit System 

We first estimate the difference between corrector (2.24)-(2.25) and the formal limit system (2.5)-(2.6). We 
will prove the following 

Lemma 2.2. For every e > 0 let F and 9^’^ correspond to solutions of (2.24)-(2.25) and of (2.5)-(2.6) 
respectively both starting from initial conditions 9g g distributed as any Rgpe, where p^ is an invariant 
measures of the system (2.1)-(2.2). Then, for any t > 0, and any e > 0, 

E\F(t)-9°g’^{t)f <C£, (2.26) 


where C = C(Ra,t, |cr|) is independent of e > 0 and (the possibly non-unique) ps. 

Proof. Set := := F — 9g’^. Comparing (2.24)-(2.25) to (2.5)-(2.6) we see that and 

satisfy 


v^(t) = ^{Ra ■ e(^^)) -|- exp(—e ^ait)u 

^+a2(r) = b2(ug^9ln - fo2(fi^r) 


ls-exp{-£ ^ait)a^^{Ra-e(9lig)), 

= -b2(v^,9Y) - b2{u^,tp ^), - 0 ^( 0 ) = 0 . 


(2.27) 

(2.28) 
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First, from (2.27) using (2.3) we obtain 

< 3i?a2|enarTlV’"WP + 3exp(-4£-4)(|«g^5p + Ra^\e\^\a-^\^\9ls\^). (2.29) 

On the other hand from (2.28) with (2.3), (2.4), we obtain < — 2 {b 2 {v^,d^g'^),fp’^) and hence 

with (2.29) we have, for any t] > 0, 


— < C\v^\^ + ^\9Y?W? < [cRa^ + 1-01^ + C'exp(-4£ 4)(|ug_sp + Ra^\9ls\^), 


where C = C(r], la^ ^|, |e|, I 62 I) is independent of e > 0. Using Grdnwall’s inequality, /J exp(e ^(s — f)) ds < e 
and 0 '^(O) = 0 we have 


< eCexp (^CRa^t + | ^ (|u§,sP + Ra^\9lg\^), 

for a constant C = C{ri, |aj"^|, |e|, I 62 I) independent of e > 0. We now take expected values and apply Holder’s 
inequality. By choosing rj > 0 sufficiently small in order to make use of (2.18) for 9^’^ and then applying 
( 2 . 21 ) and ( 2 . 22 ) to bound terms involving 9q g and Uq g we obtain 


E|'0^(t)P < eC{l +Ra^)ejip{tC{Ra^ + Icrp)) • (Eexp(?7|6»g gp) • E(|ug 
< eC(l + Ra^) exp {tC{Ra^ + |crp)) 


for a constant C = C{\ai ^|, |e|, I 62 I) which independent of e > 0. The proof is now complete. 


□ 


2.4.3 Estimates between the Corrector and Small Parameter Solutions 

We next compare solutions of the small parameter system (2.1)-(2.2) to the corrector (2.24)-(2.25) and prove 

Lemma 2.3. For any e > 0 let 9g and he solutions of (2.1)-(2.2) and of (2.24)-(2.25) respectively both 
with same initial conditions distributed as Refie; where is any invariant measure pL^ of (2.1)-(2.2). Then, 
for every t > 0, there exist C = C{Ra, \<j\,t) > 0, 5 = S{Ra, |cr|,t) > 0, both independent of e > 0 and the 
particular choice of such that 

e( sup |6'|(s)-0'^(s)pM <Ce^, (2.30) 

\se[o,t] J 

for every e > 0. 

Proof. Set = Ug — = 9g — 9^^. Similarly to (2.28), 0® satisfies 

+ a^m = e%) - 62 (^i^ r), 0"(O) = 0 , 

and hence with (2.3)-(2.4) 

+2\r? <2\{h2{i^,9%),r)\- (2.31) 

Given that '0®(O) = 0, and 1 ( 62 ( 1 )^, 0|), 0^)1 < + C|0|p|?;'^p we have that 

ptAr 

sup |0®(sAr)p<C' sup |0|(sAr)p / \v‘^{t')\^dt' (2.32) 

se[o,t] se[o,q Jo 

for any t > 0 and any stopping time r, where C = C(| 62 |). 

In order to obtain bounds for v^ we rewrite (2.24) as 

ai{vd{t)) = Ra ■ e{9‘^) + exp(—e“^ait)(ai('Ug 5 ) — Ra-e{9Q g)). 
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As such, comparing with (2.1), 


/ dul 


e + bi{us,Us)j + ai(u®) = i?a- e(V'®) - exp(-e ^ait){ai{ul g) - Ra-e{0^^ g)). 

Our next step is to derive a stochastic equation for u^. Returning to (2.24)-(2.25) we obtain 
du^ = a^^{Ra ■ e{d9^)) — £~^ exp(—[ 01 (^ 0 ,s) ~ ^“‘ 6 (^ 0 ,s)] 

N 

= — a^^{Ra ■ e{b2{u^, d‘^) + a2{d’^)))dt + ^ a^^{Ra ■ e{ak))dW'^ 

k=l 

— e~^ exp{—£~^ait) (ai(itQ g) — Ra-e{9Q g)) dt. 

Combining (2.33) and (2.34) we infer 

dv^ + -ai(v'^)dt . el’tp'^idt 

£ £ 


(2.33) 


(2.34) 


N 

(yi{ug, Ug) — a^^{Ra ■ e{b2{u^,9^) + a2(9^)))"j dt — ^ a^^{Ra ■ e{(Tk))dW^ , 


fc=i 


and hence by the Ito lemma 


N 


d\v H—(ai(u ), V )dt = - {e{'ip ), v )dt + Oi {Ra ■ e{ak))\^dt 

£ £ ^—' 

k=l 

- 2 (yi{ug,ug) - a^^{Ra ■ e{b2{u^, 9‘^) + a2{9‘^))),v‘^) dt 

N 


- 2^(ar^(Ra • e{ak)),i^)dW\ 


(2.35) 


fc=i 


Next, with (2.3) and routine bounds we infer 


1 , 


,CRa^ 


<———+ CRa'^\a\^dt + Ce^luH^ + 


dt 


N 


-2^(ar'(i?a-e(afe)),r)dlR\ 


fc=l 


where C = C(\ai\,\a^ ^|,|e|,5i) is again independent of e > 0. Integrating in time, multiplying by e and 
recalling that ^^(O) = 0 we obtain 


pt/\T 


pT/\r pt/\T ^ ^ \ 

^CRa^ J (|'0®pd<' + £\<T\‘^)dt' + Ce^ J + Ra^\u^\‘^\9‘^\‘^ + \9'^\ydt' 

^ ptAr 

-2£V/ {a^\Ra-e{ak)),i^)dW'^, (2.36) 

Jo 


for any t > 0 and any stopping time r. 

We now combine (2.32) and (2.36) as follows. Consider the stopping times 

Tk ■■= inf{|6»|(t)|^ > k}, 


(2.37) 


for each k > 0. Applying (2.32) with these stopping times, making use of (2.36) and then taking expected 
values yields 

E sup |'!/(‘^(s A r„)p </ti?a^C [ (^E sup {s A Tk)\^ + £\a\^^ dr 

sG[0,t] 


0 \ sG[ 0 ,r] 


+ £^kC / E 


e |4 


+ Ra2|u"n0"|2 + |0^|2^ds, 


(2.38) 
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where C > 0 is independent of i?a, k > 0 and £ > 0. In order the bound the final term in (2.38) we estimate 
defined in (2.24) by a triangle inequality. Next, we use that 9^ satisfies drift diffusion equation (2.16), 
and consequently we apply (2.18), (2.21)-(2.22) to infer 

E < C(1 + i?a^), 

where C = C(|fT|, |ai|, |aj"^|) is independent of £ > 0 and Ra. Hence, 

E sup |'0‘^(s A r„)p < k(1 + f /^E sup |^/;^(s A r^)!^ + 

sG[0,t] Jo \ se[0.r] / 

and hence, with the Gronwall inequality we find 

E| sup |'0‘^(s)pllr,,>i ) < E sup |^^'^(s A Tk)!^ < £t exp(KC'i) (2.39) 

\sG[o.t] J sefo.t] 

which holds for a constant Ci independent of £ > 0 and k > 0. 

Set Xe{t) := supsg[o_t] \tp^{s)\'^ and for each t>0, k> 0, £>0 define the sets 


Et,K,E ■= < sup > K > = {tk < t}. 

(sG[0.t] J 

For each t > 0 and each sufficiently small t] = 77(|cr|) > 0 one finds by Markov inequality, (2.18), and (2.21) 
that for small r] > 0 


P(£^t.«.e) < e-’^'^Eexp 


77 sup |6»|(s)| 

sG[0,t] 


< e-’"'C'exp(?7t|cTn, 


where C > 0 is independent of £ > 0 and t,K > 0. On the other hand, by (2.39) for any <5 > 0, 


/ 






XeitYH 


V 




— ^^E(Xe(t) lLi-fc<tllrfc+i>t) 


oo 

< ^ (E(Xe(t)]l.,^,>t))' (P(rfe < t))'-' 

oo 

<e^t^ exp{rit\a\^)C ^ exp((5C'i(fc + 1) — (1 — S)r]k ), 

k=0 


(2.40) 


(2.41) 


where we have used (2.39), (2.40) for the final bound. Here C is independent e and Ci is the constant 
appearing in (2.39). Thus when S < the series in (2.41) converges and the desired bound (2.30) follows 
for any such value of 6 . The proof is thus complete. □ 


3 The Boussinesq Equations and Their Mathematical Setting 

We begin our analysis of the stochastic Boussinesq Equations by recalling some details of the mathematical 
setting. The section concludes with a mathematically precise restatement of Theorem 1.1. Here and below 
we implicitly assume that C, c, Cq etc. are constants depending on the domain 2?, any other dependence is 
indicated explicitly. 

For the forthcoming analysis it is convenient to consider an equivalent, homogenous, form of the stochastic 
Boussinesq Equations. Introducing the ‘small parameter’ e = Pr~^ > 0 and making the change of variable 




Large Prandtl Number Asymptotics in Randomly Forced Turbulent Convection 


13 


9^ = T — Ra(l — z) we can rewrite (1.1)-(1.2) as 

£(atu^ + • Vu^) - Au= = + i?akr, V • = 0, (3.1) 

N 

dO^ + • VO^dt = Ra ■ u^dt + A9^dt + ^ OkdW’^, (3.2) 

k=l 

supplemented with the homogenous boundary conditions 

Up=o = = 0; ^\z=o = ^\z=i = Oj periodic in x = (xi,X 2 ). (3.3) 

Note here that we have implicitly modified the pressure in (3.1) by Ra{z — ^z^) since (1 — z)k ~\/(z— ^z^). 
The corresponding infinite Prandtl system (e = 0) is given by 

-Au° = Vp + i?ak6»°, V • u° = 0, (3.4) 

N 

de° + u° • we°dt = Ra ■ u°^dt + Ae°dt + (TkdW’^, (3.5) 

k=l 

again with initial conditions 0°(O) = 0g and boundary conditions as in (3.3). 


Remark 3.1. Notice that we do not prescribe an initial condition for u*^ as this component of (3.4)-(3.5) 
does not satisfy an independent evolution equation. Indeed (3.5) can be rewritten as 

N 

d9° + {L9°) ■ ve°dt = Ra{Le°)ddt + Ae°dt + Y <^kdW^ , (3.6) 

k=l 


where L = RaA ^Pk for A the Stokes operator and P the Leray projector, that is L9 is the solution of 
—Au = Vp + Rak9, V • u = 0; see (5.2) below. 

Remark 3.2. The systems (3.1)-(3.2) or (3.4)-(3.5) can be reformulated in terms of T = 9^ + Ra(l — z), 
which satisfies (1.1)-(1.2) or (1.4)-(1.5) and have boundary conditions given as in (1.3). Our analysis will 
make use of both of these formulations. 

Remark 3.3. As noted above, parameters in the problem are the Prandtl (Pr = e~^) and Rayleigh numbers 
(Ra, Ra), which are unit-less. In terms of basic physical quantities of interest we have that 


e ^ = Pr = —, Ra = 

K 





7 


(3.7) 


where v is the kinematic viscosity, k the thermal diffusivity, g the gravitational constant, a the coefficient 
of thermal expansion, h the distance between the confining plates, Tb — Tt the temperature differential, and 
7 = R!pc the intensity T-L of the volumetric heat flux‘d normalized by the density p and specific heat c of the 
fluid. We refer the interested reader to [FGHR] where the dimensionless form of the stochastically driven 
Boussinesq equations, (1.1)-(1.2) or equivalently (3.1)-(3.3), are derived. 


3.1 The Functional Setting 

We next define the phase space for the Boussinesq equations. Our setting is very close to the classical 
framework for the Navier-Stokes equations; see e.g. [CF88, TemOl] for further details. 

For every e > 0 define H := Hi x H 2 as the phase space for (3.1)-(3.3), where 

Hi := {u G {L^{V)Y : V • u = 0, u • ^ = 0, u is periodic in x}, 

H 2 ■= {9 G L^{V) : 9 is periodic in x} 

■^In our setting H has units of power/Vvolume*time. 
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and we abuse notation setting hy H = for (3.4)-(3.5) when £ = 0. As usual H is endowed with the 
standard L^-norm which we denote as || • ||. All other norms are written as || • |lx below for a given space X. 
We define type spaces as 

Vi := {u G {H^{'D))^ ; V • u = 0, = 0, u is periodic in x}, 

V 2 := {0 G H^{T)) : 0 | 2 =o,i = 0,0 is periodic in x}. 

Let y = Vi X V 2 if £ > 0 and 1^ = V 2 if £ = 0. Finally set 

lig ■. H ^ H 2 to be the projection onto the 6 component of H. (3-8) 

As above in (2.7), for any Borel measure p, G Pr(i7), we take Ilgfi to be the push-forward of p by Rg. 

We have the following general results concerning the existence and uniqueness of solutions of (3.1)-(3.3) 
and (3.4)-(3.5): 

Proposition 3.1 (Existence, Uniqueness and Continuous Dependence of Solutions on Data). 

(i) For every e > 0 and any given p° G Pr{H) there exists a stochastic basis S = (D, P, IF) 

upon which is defined an H-valued stochastic process with the regularity 

(u^r) G L2(f2;LL([o,oo);U)nL^,([0,oo);i7)), 

which is weakly continuous in F[, adapted to {iFt}t>o, satisfies (3.1)-(3.2) weakly and such that 
(u'^(O), 0*^(0)) is distributed as p°. We say that such a pair (5, (u®,0'^)) is a weak-martingale solu¬ 
tion of (3.1)-(3.3). 

(a) Additionally, for any £ > 0, there exists a martingale solution {S, (ug,0g)) which is stationary in time. 
These stationary solutions {S, (ug,0|)) may be chosen in such a way that, for any p> 2 there is an 

V = v{p) > 0, 

sup/ exp(? 7 (||uf-k ||0|||p))dpe(u,0) < 00 , (3.9) 

e>0 JH 

where Pe)-) = P((u|, 0|) G •)■ 

(Hi) Now consider the case when £ = 0. Fix a stochastic basis S and any iFo measurable random variable 
00 G L^{fl,FI). Then there exists a unique process 0° 

0° G LL([0, 00 ); U) n C([0, oo); H)) (3.10) 

which is J-t-adapted, weakly solves (3.4)-(3.5) and satisfies the initial condition 0°(O) = 0o. 

(iv) For a given stochastic basis S and each Oq € H denote 0°(-,0o, W) as the unique corresponding stochas¬ 
tic process satisfying (3.4)-(3.5) with (3.10). We have that 6 q 0°(t,0o, VF) is Frechet differentiable 
in 00 G H for any t > 0 and any fixed realization Wf) = W (-, 0 ;). On the other hand W 1 —>■ 0°(f, 0o, W) 
is Frechet differentiable in W from C'o([0, t],R^) to H for each fixed 0o G H and t > 0. 

These results are standard for a system like (3.1)-(3.3); see e.g. [DPZ92, FG95, DGHTll]. The only 
novel difficulty in view of existing methods is the uniform moment bound (3.9). The existence of such a 
collection of solutions is established using the maximum principle and exponential moment bounds in the 
companion work [FGHR] ; cf. Appendix A below. 

The Markovian framework for (3.4)-(3.5) is defined as above for the toy model (2.5)-(2.6). The transition 
functions are given by 

P°(0o,A) :=P(0°(t,0o) G A), t>O,0oGH,AGB{H) (3.11) 

and associated semigroup by 

P°(/)(^o) :=E</>(0°(t,0o)), t>0, cfGMbiH), (3.12) 
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where Mi,{H) is the set of bounded measurable functions on H. This semigroup acts on Borel probability 
measures p according to 

pP°{A)= [ P°{e,A)dpie), A€BiH). 

JH 

In view of the continuous dependence on initial conditions the semigroup is Feller, that is, it maps 

the set of continuous bounded functions on H, Cb{H) to itself. 

3.2 Statement of the Main Results 

We now precisely formulate the main result of the work (cf. Theorem 1.1). 

Theorem 3.1. Let {P°}t>o be the Markov semigroup associated to (3.4)-(3.5) defined in (3.12). There 
exists Nq > 0 and tjq > 0 depending only on Ra and Ra such that 

(i) if N > No, that is, if we directly force the first N eigenfunctions in (3.5), then 

Pr,id-Pt,i^Pt) < C'exp(- 7 t)p^(/x,/i), (3.13) 

for any p, p G Pri(H), rj G (0,r]o) o,nd every t > 0, where prj is defined in (2.9). In particular there 
exists a unigue ergodic invariant measure p G Pri{H) of (3.4)-(3.5). 

(ii) There exists a collection of measures {pe}e>o corresponding to statistically invariant states of (3.1)- 
(3.3) which satisfy the uniform moment condition (3.9) for some sufficiently small p > 0 and p > 3. 
For any such collection, if N > Nq we have that 

p^{PePe, do) < , (3.14) 

where C = C{Ra, Ra) and q = q{Ra,Ra) are both independent of e > 0. As such, 

nm)-^{.o%))<c\m2 ^e\ ( 3 . 15 ) 

for any element G V 2 r){H) as defined in (2.11). Here { 05 }e>o and 9g are stationary solutions of 
(3.1)-(3.3) and (3.4)-(3.5) respectively corresponding to {pe}e>o and pof 

The proofs of (i) and (ii) are carried out in Section 4 and 5 respectively. We conclude this section by 
making several important remarks. 

Remark 3.4. 

(i) Using the approach detailed here, the infinite Prandtl limit analogous to (3.14), (3.15) may also be 
established for a two dimensional version of (3.1)-(3.3). Here one can also show that (3.1)-(3.3) has 
a well defined Markov semigroup. Thus, any statistically invariant state corresponds to an invariant 
measure of the associated semigroup. This additional ingredient of Markovianity in two dimensions 
allows us to show in [FGHR] that the e-independent exponential moment bounds in (3.9) hold for all 
invariant measures. 

(ii) In 3D the existence of a sequence of statistically invariant states of (3.1)-(3.3) satisfying the uniform 
moment bound (3.9) is established in the companion work [FGHR]. Here by contrast to the 2D case 
we have not been able to show that every sequence of statistically invariant states of (3.1)-(3.3) have 
such (uniform) exponential moments. 

(Hi) An interesting outstanding issue is to establish the convergence of the velocity fields in the large Prandtl 
number limit. This convergence of the extended measures p^ as e ^ 0 as compared to (3.14) will be 
addressed in future work. 

®It is worth noting here that are only stationary Martingale solutions. As such we cannot suppose that the collection of 
these solution are all defined relative to the same stochastic basis. This subtlety will not cause us any trouble in what 

follows and we shall essential suppress this technical point in order to avoid notational confusion. 
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4 Contraction in the Wasserstein Distance for the Infinite Prandtl 
System 

In this section we establish some properties of the infinite Prandtl system, (3.4)-(3.5). These properties 
provides a sufficient condition for the contraction bound (3.13) as a consequence of a general result in 
[HMDS, Theorem 3.1], 

Proposition 4.1. There exist rjQ > 0 and Nq, depending only on Ra, Ra, such that whenever the number of 

forced modes N in (3.4)-(3.5) exceeds Nq and for any 0 < ij < rjo we have 

(a) Lyapunov structure: For all t* > 0, there exists Ci = Ci{t*,T]) such that 

E (exp(77||0°(t, 0°)f )(1 + Ilyo.tlD) < Cl exp(77(l + ARaRa)e-^/^\\9°f ), (4.1) 

for each 0q G F[ and every t G There the operator is the Frechet derivative of O^lt^Oo) with 

respect to initial condition 0^; see (B.l) and (B.8) below. 

(b) Gradient Bound for Markov semigroup: for any (p G Cl{H), and every t > 0, 0 G H 

WNPUm < Cexp(,7||0f) (^^pO(|<^(0)|2) + S{t)^P°{\\ymV)^ , (4.2) 

where i5(t) —^ 0 as t —^ oo. Here again 5 and C > 0 depend only on Ra, Ra, and tj. 

(c) Irreducibility condition: for any M > 0, 6 > 0 there is t* = t^{M,5,rf) such that for each t > t* 

inf sup r{{0,0) G H X H : pni0,9) < 6} > 0, (4.3) 

l|eo|M|eoll<M rcciSe^p°,5g^p°) 


where, as above in (2.9), C{6egP^,6g^P^) denotes the collection of all couplings of the measures Sg^Pjp 
and ^§gPt 

Proving the first condition, (a), essentially reduces to establishing a moment bound which follows from 
estimates found in [FGHR], and which we recall below in Appendix A (see Proposition A. 2). The second 
condition, (4.2), can be translated to a control problem though the use of Malliavin calculus which in our 
setting amounts to proving a relatively straightforward Foias-Prodi type estimate. Since both (a) and (b) can 
be established by methods essentially contained in previous works we relegate further details to appendices 
(see Sections A and B below). As already mentioned above, the principal novel challenge here is to prove 
the irreducibility condition (c) which we turn to next. 

4.1 Irreducibility 

In previous related works the proof of irreducibility essentially relies on the fact that the governing equations 
without the stochastic forcing have a trivial attractor which is stable under small force perturbations; see 
e.g. [EMOl, HM06, CGHV13, FGHRT13]. Here for (1.1)-(1.2) the dynamics without body forces can be 
highly non-trivial. Our approach shows that we can reduce (4.3) to a control problem though the Girsanov 
theorem and careful stopping time arguments. We believe that this strategy may be applicable to other 
problems. 

As a preliminary step we show that (4.3) may be reduced to proving a slightly simpler bound. 

Lemma 4.1. For a given TV > 0 consider (3.4)-(3.5) with N independently forced directions. If for every 
M, (5 > 0 there is a t^, = t.f{M, 5) > 0 such that 

inf P(||0°(t, 0o)|| < (5) > 0, for each t >t.f, (4.4) 

||eo||<M 

then Proposition 4-1 (c), that is, (4.3), holds for such an N and any r/ > 0. 
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Proof. For any 9o,6o G H consider the element f G C{S 0 gP°,Sg^P°) defined on cylinder sets as 

f(^ xB)= Ptieo,A) X Ptieo,B), A,Bg B{H). 

For each t > 0 and any M, 77 , 7 > 0 one has 

inf sup F{(0, 0) G H X H : prf{9,9) < 7 } 

l|eoll.l|eoll<M TeC{,Se^Pt,5g^Pt) 


> 


inf f 1(9,9) G BiX Bi : || 6 I|| + ||0|| < 7exp(-277)| 


l|eoll.l|eoll<M 


> 


(^* 0 , \9gH : || 6 »|| < min{ 7/2 • exp(-277), 1 }|) 


= „ < min{ 7 / 2 -exp(- 2 ? 7 ),l}) 

where we have used (2.8) for the first inequality. By now applying (4.4) with 5 = min{ 7/2 • exp(— 277 ), 1}) 
and the given M > 0 the desired result now follows. □ 

In order to establish (4.3) the rest of the section is therefore devoted to showing that 

Proposition 4.2. There exists an Nq = NQ{Ra, Ra) sufficiently large such that, for any N > Nq and every 
M,5 > 0, there is a t^. = t,,(M,5) > 0 such that (4.4) is satisfied. 


Proof of Proposition f.2. We first consider the modified system 
An = Vp + Ra\i9 , V • u = 0, 


N 


d9 + u - y9dt = {Ra ■ Ud + A0 — XMPN9)dt + ^ akdW’^, 0(0) = 9q , 


(4.5) 

(4.6) 


k=l 


and establish the analogue of (4.4) for 9 when N is sufficiently large.® For this purpose we consider the 
change of variable := 9 — = 9 — aW which satisfies 

dti^ + u • Vt/) = Ra ■ Ud + At/; — XjvPNi’ + (AcrW — XnPnO'W — u • V(fTlF)), tp{0) = 00- 

Taking an inner product with t/;, using that u is divergence free, the inverse Poincare inequality (see (B.14) 
below) and the bound 

||Vu|| < i?a||0-|| < i?a(||T/;|| + llaWlI) (4.7) 

which follows from (4.5) we have 

< (RaM + IIAaWlI + XM\aW\\ + ||u||||VaW|U==)||T^|| 

< C(i?ai?a(||T/;|| + ||aW||) + ||AaW|| + A^||aW|| + Ra{\m + ||aW||)||VaW||Loo)||T/;||. 

Rearranging we conclude 


dt 


{Xn - Ra{Ra + || V(tW||l=o))||t/;|| <{RaRa + Xn + i?a|| VaVF||Lcx,)|lCTW|| + || A(tW||. (4.8) 


Next, we use the fact that, with positive probability, each of ||crW||, ||VcrIF||, || AtTW|| stays close to zero over 
finite time intervals. For 'y>0,t>0, N>0 consider the sets 


■■= ■( sup ||VctVF||loo < 1 , sup ||AcrIF|| < sup ||crIF|| < 7 


se[o.i] 


se[o.t] 


se[o.t] 


1 


^{RaRa -t- Xjv F Ra) 


A 1 


®Here recall that Pjv denotes the projection onto the first N modes of the Stokes problem (with boundary conditions as in 
(3.3)) and Ajv is the corresponding largest eigenvalue in this collection. 
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Since a is spatially smooth we infer from standard properties of Brownian motion that > 0 for 

any 7 > 0, t > 0 and N > 0. On the other hand, on the differential inequality 

guv'll + (Aat - Ra{Ra + 1))||?/'|| < 7 
holds over the interval [0,t]. Hence, choosing N sufficiently large so that 

Aat > max{2i?a(i?a + 1), 1} (4.9) 


we infer that 


||0(t,0o)|| < mm + IkW^II < 2j + e-^-m\Oo\\ 

on Note that (4.9) sets the condition on Nq in the statement of Proposition 4.2. In any case, for a 

given M > 0, 5 > 0 by choosing t* = t*(M, S) such that < | we have for any t > t* 

^l^inf^P(||0(t,0o)|| < 5) > nX5,A,t,N) >0. (4.10) 

In order to now infer (4.4) from (4.10) we apply the Girsanov theorem to a slightly modified version of 
(4.5)-(4.6) in conjunction with further bounds on (4.5)-(4.6). For K > 0 and 6 q G H define 6 k = 6*if(-,0o) 
as the solution of (4.5)-(4.6) with the term —XkPnS is replaced with —XnPn6kXk{\\Pn6k\\)- Here xk is 
a smooth, non-negative cut-off function which is one for \x\ < K and zero for \x\ > K Consider the 
stopping times 


Tif( 6 »o) = inf |||PAr 0 if(t, 0 o)|| > /trj , 

for any K > 0 and any 6 q G H. It is not hard show that, for any K > 0 and any Oq G H 

P (^0(t A rK( 6 >o), 6 »o) = 0if(t A Tif( 6 >o), 6 >o), for every t > 0^ = 1. (4.11) 

On the other hand, for any do G H and K > 0, Ok{-, 60 ) is absolutely continuous with respect to the processes 
d^{-, 6 o) solving (3.4)-(3.5). Indeed for 60 G H and K > 0 define 

M9o,K(t) =exp J aeo,KdW - ^ J \aeo,K\^ds^ , (4.12) 


where 


Cteo,K = —XnCT ^ PN6K{,t,6o)XK{,\\PN6K{t,9o)\\) 


and define 


dQ9o,K '■= M9a,KdP ■ 

Notice that, since \(J~^Pn 6 k{s, 6 o)xk{\\Pn 6 k\\)\ < ’ {K + 1), the Novikov condition is satished and 

Girsanov theorem applies to 6 k{-, 60 ) and Xi 9 ^^K for any K > 0 and any 60 G H. 

Now, according to the Girsanov theorem, for any S > 0 and t > 0, we have 

P(||6»(t, 6»o)|| <S)= Q9o,K{\\6Ki.t, 6»o)|| < 5) = E {M9o,K{m\\9j^{t,9a)\\<^ ’ 

and hence for any S > 0, 60 G H and for any 2., K,t > 0, Markov inequality implies 

¥{\\6{t, 0o)|| < <5) > HP (\\6K{t, 0o)|| < <5, M 9 o,K{t) > 3 ) > 3P {m, ^o)|| < <5, M 9 „Kit) > 2 , tk{ 6 o) > t) , 
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where we used (4.11) for the final inequality. On the other hand 


P(||0(t,0o)ll <<5) <P(||0(t,0o)ll <<5,A^eo(t)>3)+P(A^«o(i) <3) 

< P (||0(t,6»o)|| < 5,Mea,K{t) > !I,rK(6»o) > t) F'P {MKfioi^) < +P(t/c(6»o) < t). 

These two bounds yield 

inf P(||0(t,0o)|| <5) 

l|eoll<M 


inf P(||0(f,0o)|| <<5)-3 sup P(7W^,e„(t) <n)+P(TK(0o) <t) 

l|eoll<M ||0 oII<m' 


(4.13) 


for any M, 6,t > 0 and for any K,2 > 0. 

Since the first term on the left hand side of (4.13) is independent of it' > 0 and has the same dependence 
on H > 0 as the second term we finish the argument by showing that for every fixed M, K,t > 0 


and for every given M, t > 0 


sup P (A4eo,K(f) < 3) —>■ 0, as H —>■ 0, 

l|eoll<M 


sup F{TK{do) < t) —>■ 0, as itT —>• oo. 
11^0 ll<M 


(4.14) 


(4.15) 


For the first bound (4.14), we have from (4.12) and Ito isometry 

P {Mea,Kit) < H) =P aeo,KdW J Ic^So.Kl'^ds > log(!l~^) 

+ ^y \aeo,K\'^ds 


< 


< 


log(3-i) 

2 


E 


oteo,KdW 


log( 3 _i)lE(^l + A^||a-^P^ry^ \\PN0{t,0o)rxK{\\PNOK{t,9o)\\)ds 
2{l + X%\\a-^PNf{K + l)H) 


log(:i-i) 

valid for any 3 G (0,1), K > 0, and any 6q G H. For the second bound observe that, in view of (4.11), 


^(^^(^o) < 0 < P 1 sup ||PAr0(s,6>o)|| > if ) < ( sup \\0{s,9o 


se[o.i] 


se[o.t] 


(4.16) 


From the Ito formula follows 

d\\e\\ + 2XN\\PN^^dt + 2\\vefdt = {2Ra{ud, 9) + ij dt + {a, OjdW. 

Integrating in time and using (4.5), inverse Poincare inequality (see (B.I4)), and (4.7) we infer for any s > 0 


||6»f + 2AAfE/ ||6»fdr < ||6lof+ 2i?ai?aE / ||6lf dr+ s +2 sup 
Jo Jo r-e[0,s] 

Using the assumption (4.9) and the Birkholder-Davis-Gundy inequality we infer 


(cr, 9)dW 


E sup ||6»(s,6»o)f <\\9of + m. 

Vse[o,t] / 


(4.17) 


Combining (4.16) and (4.17) thus yields the second bound (4.15). 
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With (4.10), (4.13), (4.14), and (4.15) now in hand we now conclude the proof by arguing as follows. 
Given any <5 > 0 and any M > 0 choose as above so that < |. Fix any t > and by 

(4.10) a = a{M,6,t) := inf||ep||<MIP(||0(i,^o) < d) > 0. Now by (4.15) we can pick K sufficiently large 
so that suppg||<jvf P(T/f (0o) < t) < a/4. With K,M,t fixed we choose H > 0 small enough so that 
sup||gjj||<jvf P (.Adi<'_6(o(^) < ^ a/4- Finally by combining these choices with (4.13) we hnally obtain that 

inf ^P(l|6'(^:fi'o)ll < <5) > ^ > 0. 

||6lo||<M 2 

The proof of Proposition 4.2 is thus complete. □ 


5 Convergence on Finite Time Intervals 

Having established the contraction condition (3.13), we now prove the second part of Theorem 3.1. Observe 
that (3.13) implies an analogue of (2.13) for (3.4)-(3.5), where | • | is replaced by the norms || • ||. In view 
of the uniform bound (3.9), Theorem 3.1, (ii) is thus an immediate consequence of the following proposition. 

Proposition 5.1. For each e > 0 let be the second component of a statistically stationary solution of 
(3.1)-(3.2) satisfying the uniform bound (3.9) and let be solution of (3.4)-(3.5) with the initial condition 
^S^(O) ■= ^so- Then, for any e > 0 and any t>0, 

E||0|(t)-0^’"(i)f (5-1) 

for constants C = C{Ra,Ra,t) > 0, i5 = S{Ra, Ra,t) > 0 independent of e > 0.^ 

5.1 Defining the corrector 

In order to proceed as in Section 2.4.1 and define the ‘corrector’ we first recall some properties of the steady 
Stokes problem 


-Au = Vp + f, V-u = 0 (5.2) 

and the associated linear evolution given as 

e9tu — An = Vp + f, V • u = 0, u(0) = uq, (5.3) 

for any e > 0 and relative to the (sufficiently regular) data f, uq. Both (5.2) and (5.3) are supplemented 
with same mixed periodic-Dirichlet boundary conditions as for u in (3.3). 

Classically, the equations (5.2) and (5.3) can be understood in terms of the functional framework given in 
Section 3.1; see e.g. [TemOl] for a systematic treatment. Recall that the Stokes operator is given by A = —PA 
which is self-adjoint, positive and unbounded on Hi. Here P is the Leray projection on divergence free vector 
fields {Lf{T>)Y Hif As such we may rewrite (5.2), as Au = Pf. Along the lines of the classical elliptic 
theory one may show that for any f G (L^(T>))^ there exists a unique u € P(A) = Vi fl denoted 

by u = A~^Pf which satisfies 

\\A-^Pf\\H^<C\\{\\. (5.4) 

Turning to (5.3) we observe that for any f € oo); Hi) and uq € Hi there exists a unique solution u 

of (5.3) relative to this data with u G ^/^^([O, oo); Z1(A)) fl (^([O, oo); Hi). The flow associated with (5.3), i.e. 
(5.3) with f = 0, defines an analytic semigroup which we denote as exp(—e“^Af). 

^Indeed, given 9g and its associated basis S'^, we can obtain (a unique) pathwise solution of 6^’® relative to this basis 
according to Proposition 3.1, (ii). While the forthcoming estimates may therefore take place in a different stochastic basis for 
different values of e > 0 the constants C and <5 will be shown to be independent of particular sequence of basses. 

^Equivalently Au = —Au — Vp, where p = p(u) the ‘pressure’ is the unique function satisfying Ap = div(Au) in the 
weak sense. 
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With these preliminaries in hand, essentially following (2.24)-(2.25), we define the corrector system. 

iF(t) = exp{—e~^At)P]veUg Q + A~^{P{Ra ■ k0'^)) — exp{—e~^At)A~^{P{Rak0g Q)) 

= A-i(P(i?a-k0'^)) + w^(t), (5.5) 

N 

dr + (u'^ • - A0") dt = Pa • ^ r(0) = 6»|^o> (5-6) 

fe=i 

where w'^(t) = exp(—e“^At)(Pjveu| q — A~^{P{Ra ■ kd| q))) so that under this definition u®(0) = PatsuI g. 
Here, for technical reasons, we slightly modify the initial condition on u® in (5.5) compared to (2.24) by 
taking such that ~ 1, where P/ve is the projection onto the first N’^ modes of the Stokes operator 
A. Observe that for t > 0, u® solves 

-Au'' = Vp^+Pa-kr+ Aw^(t), V-11^=0, (5.7) 


with solving 

edtw^ — Aw^ = Vg®, V • w® = 0, 

(5.8) 

w^(0) = PA/eul o - y| o, where y|_o solves Ay|_o + Vg = Pa • k6»|_o, V • y|_o = 0 . 
and supplemented with the same boundary conditions as u in (3.3). 

5.2 Comparing the corrector to the infinite Prandtl System 

Having defined the corrector (u®, d^) we proceed as above to prove (5.1) in two steps. We begin by estimating 
difference between the corrector and the infinite Prandtl system. In comparison to Lemma 2.2 we have 

Lemma 5.1. For e > 0, let (u®,d®) be the solution of (5.5)-(5.6) and let (ug^,dg^) to he the solution 
of (3.4)-(3.5) with initial condition 9gQ being the second component of stationary solutions of (3.1)-(3.2) 
satisfying (3.9). Then for any e > 0 and any t > 0 there exist C = C{Ra, Ra,t) and 6 = 6{Ra, Ra,t) both 
independent of e > 0 such that 

]E||r(t)-d°’^(t)f (5.9) 

Proof. Let and cjf = 9^ — Using that 9^ and dg^ share the same initial condition we see 

that and (j)^ satisfy 

-Av''= + Pafc((i^ - Aw^ V-v''=0, (5.10) 

- A(j)^ =Ra-vl--v^-^9°^^-iTf{Q) = 0, (5.11) 

where w'^ obeys (5.8). 

Starting from (5.11), using that W, v® are divergence free along with the Poincare inequality and standard 
Sobolev embeddings we estimate 

dx < P'a||<^^||||v-|| + ||v^|Ua||V</)^||||d°’^|h3 
< + CllVv^f (pV + ||d°>^||23) . (5.12) 

Next, testing (5.10) with v'^ and using Holder and Poincare inequalities, we find 

||Vv"f < 2Ra^U^\\^ + 2||Vw"f. (5.13) 

To bound we infer directly from (5.8) that 

|l|w^(t)f + ||Vw^(s)f ds < < |(||u|,g|| + Pa||d|,g||)^ 


(5.14) 
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for every t > 0. 

Rearranging in (5.12), applying (5.13), (5.14), and recalling that (/>®(0) = 0 we infer 


sup ||(^''(s)|P < C sup + ||6's’^(s)||i3 

se[o.t] se[o.i] ^ 


Ra^ sup ||(/)®(t')frfT’ + e(||u|._o|l+i?a||0|_o||)" 

t'G[0,r] 


and hence with GronwalPs inequality we have 


sup ||?i>^(s)f < e(||u|^oll + -Ra|l^l,oll)^ exp 

sG[0,t] 

Choosing 6 = S{Ra, Ra,t) sufficiently small we infer the desired bound, (5.9) by applying Proposition A.l 
and using the condition (3.9). □ 

5.3 Comparing the Corrector to the Large Prandtl System 

The harder case is to compare the corrector to the large Prandtl system. The reader may refer back to 
Lemma 2.3 for the for the analogous bounds for the toy model. 

Lemma 5.2. For each e > 0 suppose ® stationary solutions of (3.1)-(3.3) satisfying the uniform 

bound condition (3.9). Take to be the solution of (5.5)-(5.6) with the initial condition 0% q = 0s(O). 

Then for any e > 0 and any t > 0 there exist C = C{Ra, Ra,t) and 5 = S{Ra, Ra,t) both independent of 
£ > 0 such that 

E||r(<) - 6»°’^(t)f < Ce^. (5.15) 

Remark 5.1. Note that C{Ra, Ra,t) —t oo as t ^ 0+. This is due to the fact that unlike u^, u® does not 
satisfy an evolution equation and we do not prescribe an initial condition for u® in (5.5). 

Proof. Define = u| — u® and ff = 0^g — 0^. Referring to (5.5)-(5.6) and to (3.1)-(3.3) we see that ff 
satisfies 

dt^‘^ — = Ra ■ — ■ V0^ — Ug • V^'^, ^“^(0) = 0 , 

and therefore testing with and using that V • = 0 we have 

+ II= / {Ra • - V . V0^Wdx < i?a||v^||||<^^|| + ||v^|he||V^^||||r|h3 . 

Hence from standard Sobolev embeddings and the Poincare inequality 

J^\m‘^<c{\m\h+Ra^) liVv-f. 


C(i?a^t + 1) 


sup 


'Ra^ + \\e°/{s 


li.)). 


Integrating in time we infer that 


ptAr 

) \\yr{t’)fdt’ (5.16) 

for any t > 0 and any stopping time r > 0. 

We now turn to derive an evolution equation for v^. Recalling that u'^ satisfies (5.7) and comparing this 
equation with (3.1) satisfied by u| we find (cf. (2.33)) 


sup ||^^(s A T)f < sup (||r(s)||^3+i?a^ 

sG[0,t] s€[0,iAT] ^ 


£{dtUg + u| • Vug) — Av® = Vg® + Rakcp'^ — Aw^, 


( 5 . 17 ) 


Large Prandtl Number Asymptotics in Randomly Forced Turbulent Convection 


23 


where, as above, maintains (5.8). From (5.5), (5.6), and (5.8), we find that satisfies 
= 5tw^+ RaA-\P(kd0^)) 

= -(-Aw^ + Vg'^)-i?aA-i(p(k(u" • V^- - Ra ■ + RaJ2^~\Pi'^^k))dW^ . (5.18) 

^ fc=i 

Multiplying (5.17) by e~^, subtracting the resulting system from (5.18) and rearranging we obtain 


dv^—-Av^dt = - (yq^ + Rak^^'j 


dt 


N 

+ Ra (^A-^(^p(k(u^ ■ Vd" - Ae^ + Ra ■ uijjj - u| • Vu|)dt - RaJ2^~^(Pi^<^k))dW'^ , (5.19) 


k=l 

with V • = 0. Using (5.19) we now estimate as follows. The Itb formula and (5.20) reveals that 

2„__,„.. 2 
e 


=-Ra{(^‘^,v^)dt + 2Ra (a • Vd'^ — AO^ — Ra ■ '“d))) “ Ug • Vug, v® ) dt 


N 


N 


+ Ra^ J2 \A-\P{kcrk))\^dt - 2Pa^(A-i(P(kCTfe)), v^)dVF'= 

fc=i fc=i 

:=(Ti + r2 + Ta + r4 + Tg + Te)dt + SdW. 

With the Young and Poincare inequalities we have 

For T 2 we use that A~^ is self-adjoint on id, D{A) C H and that u®,v^ are divergence free, to obtain 


(5.20) 


(5.21) 


IT 2 I =2Pa 


• Ve^{A-^v^)adx 


= 2Ra 


■V{A-^v^)d9^dx 


where {A represents the third component of the vector field A Hence (5.4) and the imbedding 

^ imply 


Vv-ll < -liVv-f +eC(Pa4||rr + llfi^f). 


IT 2 I < 2Pa||fi^||||r||||V(dl-iv-)|U^ < CPa||fi^||||d~ 

Using (5.5), the triangle inequality, and (5.14) we have 

lltill < CRa\m + l|w^(i)|| < ||u|,o|| + CRamj + ||r||). 

Combining this observation with the previous bound we infer 

IT2I < ^llVv^f + sCiRa'^imf + ll^l.oll") + llul.oll") • 

For the terms T 3 and T 4 we use the regularity of Stokes operator and (5.22) to obtain 

IU 3 I + |r4| < ^llVv^f + 4eRa^{Wr + ) 

“ + 1)C'(||0'^|P + 11^1,olP + 11^1,olP) ■ 

To address Tg we take advantage of an additional cancellation. Since Ug = -f we find 

|Tg| =2Pa|(u|•Vu^v=)| =2Pa|(u|•Vv^u=)| < ^||Vv®f-f 4ePa2||u®|||oo||u|f . 


(5.22) 


(5.23) 


(5.24) 
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Next, using (5.5), (5.4), the imbedding ^ L°°, and standard properties of analytic semigroups we obtain 

llu^llioo < CRaWe^^W + ||w®(t)||L=o < CRa\\9^\ + ||PAr,u| o ||^2 + ||y|^ollff^ 

<CRa\m + \\PM^uy\\H^+Ra\\ey\\ 

Recalling that is chosen such that eXj^e ^ 1 we conclude with the generalized Poincare inequality and 
further standard manipulations that 


\n\ < ^llVv-r + CsRa^iRa^Wr + ll^l.olr) + ll^^.u|,o||^.)l|u|||^ 

< ^llVv-f + CRa\m\^ + ||0|,or + llul.oll" + l|u|r). (5.25) 

Finally we observe [Tel < CRa?. Combining the bounds (5.21)-(5.25) and rearranging in (5.20) we find 

d\\rf + -jyv^^dt < + Ra^\\ar + C{l + Ra^){l + Ra^){m^+ ||0|,or+ ||u|_or+ ||u^5r+ 1) 

N 


-2RaY,{A-\P{kak)),v^)dW\ 


k=l 


where the constant C > 0 is independent of Ra,Ra and e > 0. By (5.5) we have that v'^(O) = Qn^u^q, 
where Qate := I — P^s and Pve is the projection onto the first iV® eigenfunctions of the Stokes operator. 
Consequently for any f > 0, any stopping time r 


ptAr 


ptAr 


Vv^ll dt <£:|lQAreu|_Q|F + 4i?a' 


iu^r+s\\ar)dt' 

ptAr ^ 


+ eC'(Ra'‘ + l)(Ra^ + l) 




^ ptAr 

-eRaJ 2 {A-\P{kak))A^)dW\ 


fc=l ' 


iuir 


iui,or+i 


dt' 


(5.26) 


where C is independent of e > 0, i?a, Ra. 

Next for any k > 0 define the stopping times 


Tn := inf {||d''(t)||i3 > . 


(5.27) 


From this definition, (5.16), (5.26) we now infer 
E sup ||0®(s A t)|P 

sG[0,t] 


<4Ra^ {n + Ra^) I E sup ||(;i^(s A r)f dt'+ e (k + (E||u| pf + i?a^||af t) 

' lo sG[0,t'] ^ 7 


+ eC(« + Ra^ + l)(Ra^ + l) / E + rf + 1 


dt, 


which implies with the Gronwall inequality that 


E sup ||^®(s A t)||^ < eexp ("c'(i?a"‘+ 1 ) (k + t + l') Al£(t) 
sG[0,t] ^ V / / 


(5.28) 


Me{t) := E||u|o|p 


f 


E 


+ l|ublr +11^11 +11^1,of+ lkll 


ds, 


where 
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and the constant C is independent of k, e, Ra, and Ra. 

Now in view of (3.9) and making another usage of Proposition A.l, we observe that Me is bounded 
independently of e > 0. Therefore (5.28) is of the form (2.39). As above we define 

Et,K,e ■= < sup ll^'^llia > K i = {Tk < 0 

(sG[0,t] J 

and as in (2.40) we infer from (A.3) that 


nEt,.,e) < Ce-"" (5.29) 

for a constant C = C{Ra,t) independent of £ > 0, k > 0. Thus, applying (5.29) to (5.28) as in (2.41) we 
now infer (5.15) completing the proof of Lemma 5.2. □ 

A Appendix: Moment Bounds For Stochastic Drift-Diffusion Equa¬ 
tions 

In this appendix we collect some moment bounds proved in [FGHR] which have been used extensively in the 
analysis above. 

As in [FGHR] we consider the following class of stochastic divergence-free drift diffusion systems 

N 

d^ + v ■ = {Ra ■ V3 + A^)dt -I- akdW^, C(0) = (A.l) 

k=l 

evolving on the three dimensional domain V = [0, x [0,1]. Here Ra > 0 is a fixed parameter and 
V = {vi,V 2 ,v^) any sufficiently regular and adapted, divergence free vector field. Both v and ^ are supposed 
to satisfies the same mixed Dirichlet-Periodic boundary condition as in (3.3). Recall that by the 

change of variable T = .^ -I- Ra{l — z) we may reformulate (A.l) as 

N 

dT + v ■VTdt = NTdt + ^OkdW'^, T{Q) = Tq = + Ra{l - z), (A.2) 

k=l 

where v and T satisfy boundary conditions as in (1.3). As such, bounds for ^ solving (A.l) immediately 
translate to bounds for T. 

In [FGHR] we prove: 

Proposition A.l. Suppose thatv & a.s. and is Et-adapted. Fix 

any p >2 and any initial condition € H H L^(T>) which is Eq measurable with 

Eexp(77l[^o||ip) < oo, 

for some rj > 0. Then there exists r]Q = 770 ( 0 ", Ra,p) > 0 such that for any t > 0 and any positive rj < 770 , 

Eexp (^ ^pJ’ 2+2 <CiEexp(j]\\^o\\lp+'qpt{\\a\\lp+2P/^{4Ra^ + l))j (A.3) 

for a constant C = C{Ra,p) independent oft, rj, ^q, and v. Furthermore, 

]Eexp ’ A-4) 

where again C = C{Ra,p, j]o"jjLP,21) and k = K{Ra,'D) > 0 are independent oft, rj, ^q, and v. 

We now return to the infinite Prandtl system (3.4)-(3.5) and recall a bound analogous to (A.4) but which 
uses more of the specific structure for the advecting velocity field. 
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Proposition A.2. Fix an initial condition 9 q € FI which is iFo measurable, and let 0^ = denote 

the corresponding solution to (3.4)-(3.5). There is a universal constant rj* > 0 such that for any t > 0 and 
r] G (0, ? 7 *], there exists C = C{Ra,Ra) > 0 such that 

E(^exp(^r7||6»°f+ < Cexp (77(1 + 4 i?ai?a)e-‘/ 2 || 6 l°f) . 

The proof of Proposition A.2 can be found in [FGHR]. 

Remark A.l. Using Proposition A. 2 and (B.17) below we can easily establish the Lyapunov bound (4.1) 
with 


Cl = exp 


LCRa^e^"/'^ 

V 0^ 


RaRa ) . 


B Gradient Estimates On the Markov Semigroup 

In this section we establish the gradient bound for the Markov semigroup generated by (3.4)-(3.5) in order 
to prove (4.2). For this purpose we begin by briefly recalling how (4.2) is translated to a control problem 
through the use of Malliavin calculus. We refer to e.g. [Nua09] or [NP12] for further general background on 
this subject and to [HM06, HMll, FGHRT13] for the application of this formalism in a setting close to ours. 
Define the random operators 


Jo,t^ 


lim 0°{t,0o + S^,W)-0°it,0o,W) 
(5->0 S 


(B.l) 


for any ( € H and 


Ao tw := lim 
’ 5->.0 


e%t,eii,wF 5- e^{t,eii,w) 

5 


(B.2) 


for any w G L^(D; L^([0, t]; R'^)). Here Ao,tU> = {'D9°,w), where the unbounded operator 2) : L‘^{Lt;H) >->• 
L^{n; L^{0, t, M^) ® F[) is the Malliavin derivative and w is any element in the domain of the dual operator 
(5 of 2). 

For our purposes it is sufficient to recall that any J't-adapted process in G L^(H; L^([0, t]; R'^)) belongs 
to the domain of S and S{w) corresponds to the Ito integral of w so that 


E{'SX,w) =E 



(B.3) 


for any X G Dom(2)) and any J^j-adapted w. This is a special case of the Malliavin integration by parts 
formula. We furthermore recall that 2) satisfies a chain rule namely that if </) G C^{H) and 9 G Dom(2)) 
then 4)(9) G Dom(2)) and 


D(j){9) = X(j){9)D9. (B.4) 

Gombining (B.3)-(B.4) and making use of the Ito isometry we infer that, 

XPU{Oo)C =E (V(/)(0°(t, 0o))JoA) = E 0o)) wdw^ + E (V(/)(0°(t, 0o)) {JoA - •A>,tw)) 

^ + ^/pnw^mv) (EiiJo.t^- Ao,twrf" (b.s) 

for any (j) G Cl{H), 9o € H and any (adapted) w G L^(r2; L^([0, t];R^)). 








Large Prandtl Number Asymptotics in Randomly Forced Turbulent Convection 


27 


Our desired bound (4.2) follows from (B.5) if, for every ^ G H with ||^|| = 1 there is (adapted) w = 
w{^) S L^([0,oo);R'^) such that 

- AoMOW^ < Cexp{2r^\\eo\\^)6{t), (B.6) 

nOC 

sup E / \w{^)\'^dt < C'exp(27y||0o|P)- (B.7) 

ll«ll=i Jo 


where (5(t) —>■ 0 as t —>■ oo and C, r], and 6 are independent of Bq. 

To solve the control problem (B.6)-(B.7) we observe that (B.l) and (B.2) admit explicit characterizations 
as linearizations of (3.4)-(3.5). For any f G we let p{t) = p{t,^) := which satisfies 

dtp + u° • Vp + v° • V6»° = Ra-v^ + Ap, -Av° = Vp + Rakp, V • v° = 0, p(0) = C, (B.8) 

supplemented by boundary conditions as in (3.3).® On the other hand, setting p := Ao^tw for any w G 
L^([0,t],M'^) we have 

N 

dtp + u® • Vp + V® • V0® = Ra ■ u® + Ap + ^ OkWk, —Av® = Vp + Rakp, V • v® = 0, p(0) = 0, (B.9) 

fc=i 


again with boundary conditions as in (3.3). 

Denote p{t) = p{t, ^,w) = p — p and v := v — v for any w G T®([0, oo); and ^ G H. We now choose 
w as a function of ^ as follows. Let Pn be the projection on the first N eigenfunctions of the Laplacian 
with boundary conditions as in (3.3). Set w{t) := a~^XPNP, where A > 0 and N will be selected below. 
Relative to this choice of w = w(^), p satisfies 


^tp + u° • Vp + v° • V6»° = Ra • i;)] + Ap - APatP, Av®=Vp + Rakp, V • v° = 0, p(0) = ^. (B.IO) 


Testing (B.IO) with p and v® respectively, and using that both u® and v® are divergence free vector fields, 
we obtain 

^llpf + 2||Vpf + 2A||P^pf = 2^ (Rau® - v® • V0®) pdx (B.ll) 

and 


||Vv®||<Ra||p||. 


With standard Sobolev embeddings and (B.12) we have, for any p > 0, 




<l|v°|ke||V0®||||p|U3+Ra||v®||||p|| 
<C||Vv®||||V0®||||pr/®||Vpf/2 + i?a||Vv®||||p|| 

<CRa\\V0°\\\\p\n^p\\^A + EaRa\\pf 

<||Vpf + (C(Ra)4/®||V0®r/® + RaRa)||pf 
<||Vpf+ (p||V0®f+ C')||pf, 


(B.12) 


(B.13) 


where C = C(Ra, Ra,r]) = —h RaRa and (7 is a universal constant. Also since Rv and —A commute 

we have for Qjv := I — Pn 

llVpf = -{PnP, APnp) - {QnP, AQnp) = || VR^pf + || VQ^pf > \\VQNpf > A^||Q^pf , (B.14) 

® Notice that (B.8) can also be written as 

dtp + (L6»0) ■ Vp + (Lp) - V6)° = Ra{Lp) + Ap, p(0) = C- 

where L = RaA~^Pk and A is the stokes operator, P; cf. (5.2) and (3.6) above. Similar formulations can also be given for 
(B.9), (B.IO). 

^^Of course the choice of N will determine the number of modes subject to stochastic perturbation. Observe that w is well 
defined as is the set of the first N (nonzero) eigenvectors of the Laplacian. 
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where the last inequality follows from the generalized Poincare inequality. Choose 2X = Xn (with N to be 
chosen below) and combine (B.ll) and (B.13) to infer 

j^\\pr + {XM-MyoY + c))\\pr<o, 


and hence, since p{0) = 


||p(<)f < IlCf exp ^ 770 ^ \\Ve°fdr+{C-XN)t^ ■ (B.15) 

Applying Proposition A.l we conclude that, for any 0q G H, and ry G (0,7yo], 

]E||p(t)|P < CllCf exp (? 7 || 6 l°f + {C + r]- XN)t) , 

where C = C(Ra, i?a) is independent of ^ and 0° and t > 0. By now choosing N large enough such that 
Xn > 2(C + Jyllcrlp) we obtain 


]E||p(t)f < C||^f exp , (B.16) 

where C = C{Ra,Ra) is independent of ^ and 0g and f > 0. This yields the first bound (B. 6 ). 

To obtain the second desired bound, (B.7), we use (B.16) to estimate 

nOO nOC 

E / =||a-if A^E / WPnpV dt < Cexp(77||0°f), 

JQ 7o 

where C = C{XN,Ra,Ra) is independent of 0g yielding (B.7). The bound (4.2) now follows. 

Remark B.l. ITe can use the same argument leading to (B.15) to show that 

||p(t)f< lief expiry ^‘||V 0 °fdr + Ct) . 

That is, for any rj > 0, 

lUofl <exp(^ 7 y^ \\V9°fdr + Ct^ , (B.17) 

where, as above, C = C{Ra,Ra) = —\- RaRa. 
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